A mapping approach to synchronization in the "Zajfman trap": stability conditions 

and the synchronization mechanism 
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We present a two particle model to explain the mechanism that stabilizes a bunch of positively 
charged ions in an "ion trap resonator" [Pedersen etal, Phys. Rev. Lett. 87 (2001) 055001]. The 
model decomposes the motion of the two ions into two mappings for the free motion in different 
parts of the trap and one for a compressing momentum kick. The ions' interaction is modelled by 
a time delay, which then changes the balance between adjacent momentum kicks. Through these 
mappings we identify the microscopic process that is responsible for synchronization and give the 
conditions for that regime. 

PACS numbers: 39.10.+j, 45.50.-j 



I. INTRODUCTION 

Ion traps are widely used to store ions for long periods 
to obtain high resolution spectroscopy. Ions of the same 
charge repel each other through their long range Coulomb 
potential, so the trap has the purpose of localizing the 
ions in the measurement region. Recently a quite sur- 
prising behavior was discovered [Q, where, under special 
conditions, ions in an ion trap resonator do not diffuse 
into the whole trap but stay together as a bunch for arbi- 
trarily long times. This synchronized, collective motion 
occurs only for certain parameters of the trap fields, but 
in these regions it is stable. 

This observation, which is puzzling and of intriguing 
scientific interest in its own right, has important techno- 
logical applications. Pedersen etal have suggested the 
use of the small " table top" ion trap as a time-of-flight 
mass spectrometer. As the observation time, and there- 
fore the effective length of the spectrometer, are in prin- 
ciple only limited by collisions between the ions, mass 
resolutions can be envisioned which are otherwise only 
achieved in storage rings pL Indeed, the high resolution 
spectroscopy suggested in 1 1[ has now been achieved Q . 

Until now this synchronization effect has not been fully 
understood; from the experimental observations |^ it 
appears that the trap has to be operated in a regime 
where the ions' periods in the trap increase with their 
energy. Other tests point out that the focussing of the 
beam inside the mirrors is important. Numerical simula- 
tions confirm these empirical findings, but they, too, can 
not decide, if the requirements found so far, are complete 
and if they really stem from the underlying microscopic 
process. 

A macroscopic explanation of the synchronization ef- 
fect in terms of a " negative mass instability" was recently 
presented in It confirms that the ions' period has 
to increase with their energy and shows that a minimal 
density inside the bunch is necessary to support synchro- 
nization. But, being a mean field treatment, it can not 
give a detailed explanation of the underlying microscopic 
dynamics and it is insensitive to certain properties of the 
trapping field. 



The aim of this paper is therefore, to set up a comple- 
mentary microscopic model of the ions in the trap, which 
is simple enough to be understood completely. With this 
model we then can explain the basic mechanism, deter- 
mine the necessary conditions for synchronization, check 
their completeness and finally understand how the size 
of the bunch depends on the parameters of the trap and 
those of the ion beam. 

In this paper we concentrate on deriving the conditions 
for stability and we explain and illustrate the fundamen- 
tal mechanism. The stability limits and the behavior of 
the macroscopic bunch will be presented in a following 
paper [§. 

This paper is organized as follows: I n se ction ^ we ex- 
plain the model of the trap. In section III we lay out the 
framework of the dynamics in terms of the mappings for 
two non interacting ions. Then, in section IV, we explain 
how their interaction is incorporated in the mappings. 
With the interaction added we derive the conditions for 
bunching, sectio n M , and expla in th e underlying mecha- 
nism in section |VI| . In section VII we confirm that the 
conditions and explanations given are in fact applicable 
to the experiment. 



II. THE MODEL OF THE TRAP 

The experimental setup and the observed behavior is 
described in great detail in We will therefore only 

cite what we need to build up the model. In the exper- 
iment a bunch of ions is injected into the "ion trap res- 
onator" and its width observed, when it passes through a 
ring shaped pickup electrode in the trap's center. When 
the field gradient in the electrostatic mirrors is below a 
certain threshold, the bunch does not diffuse. 

Numerical simulations with many ions performed by 
Pedersen etal have shown synchronization in a one di- 
mensional model . We will consider a one dimensional 
model, too, and use the same simplified trap potential, 
which consists of the central field free region of length 
L and two slopes with a constant gradient F, the elec- 
trostatic mirrors. The two einzel lenses, which focus the 
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FIG. 1: Model of the trap: the simphfied potential V{x) (|l|) 
with the ions taking time Tf through the central part and Tm 
to return from the mirror. The ring shaped pickup electrode 
is located in the center of the trap. 



beam into the mirrors are neglected. They are only 
necessary for keeping the beam inside the trap. This 
setup is modelled by the potential 
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We will not treat the whole bunch, but will look at 
the behavior of only two identical ions in the trap and 
explain how these two ions synchronize their motion. 

Since synchronization is a property of the relative co- 
ordinate, we describe these two ions in relative and center 
of mass (CM) coordinates: 
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Capital letters denote CM properties and lowercase let- 
ters are used for the relative coordinate. With the ion-ion 
interaction W{x) the Hamiltonian now reads: 

H=^ + ^ + qViR+^) + qV{R-§) + q^Wix) (3) 

Without loss of generality we set the charge of the ions to 
q = +1 in the following. Though W{x) could be an arbi- 
trary two body interaction we will think of it as the repul- 
sive Coulomb potential between two positively charged 
ions. 

Our model will not be exact. In order to check the 
validity of the various approximations which we per- 
form, we insert the values from the original experiment, 
in which the bunching effect was observed initially [Q. 
These values are given in table |. If not otherwise stated 
we use atomic units in the following. 

An ion with the momentum pq = \/2mEQ will spend 
time Tjn in the mirror, which can also be expressed in 
terms of the CM momentum: 

F F 



(4) 



Ions with a velocity pQ/m = P/M need the time Tf to 
pass through the central field free region of the trap: 

Lm LM 
Tf = := aTrr. 



PQ 



P 



(5) 



quantity 


value in a.u. 


length L 
mirror field gradient qF 
ion mass m (Ar"*") 

ion energy Eo 
ion momentum pi 
beam radius d 


200 mm = 3.78 x 10^ a.u. 

gQ fceV ^ ^ 5g 

m a.u 

40 amu = 7.35 x lO" a.u. 
4.2 keV = 154 a.u. 

4.76 X 10^ a.u. 
10 ^lm = 1.89 X 10^ a.u. 


time in the mirror Tm 
time in the center Tf 
dispersion parameter a 


1.48 ^s = 6.10 X 10"' a.u. 
1.41 ^s = 5.84 X 10"' a.u. 
0.956 



TABLE I: Parameters and sizes used in our calculations. 
These values are the same as in the original experiment, in 
which bunching was originally discovered 



Here we have introduced the "dispersion parameter" a = 
'"^.f as the ratio of the two times Tf and Tm- With these 
the total time T for the period of one ion becomes: 

2Lm 



T = 2T„, + 2Tf = 2(1 + a)T, 
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From this equation we calculate the "dispersion" of 
the trap: 



dT 
dpo 
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In the experiment and in the theoretical model of ref- 
erence the ions were found to synchronize, when the 
dispersion is positive, i.e, when an ion with a higher mo- 
mentum (energy) takes longer to complete one period in 
the trap. With (^) we have a < 1: the ions then spend 
more time in the mirror than in the central region. 

Two ions with different momenta have different peri- 
ods, which in turn can be reformulated as a difference 
Aa; in their distance when we propagate both ions for 
the same time T. To calculate Ax, we linearize T (||) 
around the mean momentum With the difference in 
the ions' momenta of Ap = 2p (g) the difference between 
their periods is AT = ^^Ap. During that interval the 

ions move with the velocity leading to an increase in 
separation of 



Aa;(p) 




(8) 



Here we have used (^ — (0). Again a appears as the 
central quantity, determining the sign of Ax. 

In ||l| the energy spread of the ion source is given as 
< 0.1%, which translates into a relative momentum 
of p < 4 a.u. With the parameters of table | the rel- 
ative error between the linearized AT and the correct 
difference T(f + p) - T(f - p) is well below 0(10-^). 
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Pedersen etal defined a parameter a, too, which we 
refer to as ap [p|. From its definition as ap = we 
calculate with 



(x',p'< p) 



the following relation: 

I — a 
2Eoap = — — 



(9) 



It should be noted that yields essentially the same 
information as about the sign of the trap's dispersion. 



III. TWO NON INTERACTING IONS 

Before we look at the coupled motion of two interact- 
ing ions, we will study the non interacting case, where we 
set W = in (^) . Adding the interaction will then mod- 
ify this model and the relevant mechanism will become 
clearer. 

Instead of solving the actual equations of motion for 
the two ions we "take apart" the trap and follow the evo- 
lution of the relative coordinate through the three regions 
of the trap (cf. figure |l|): (i) the central field free part, 
(ii) the mirrors and (iii) the kink region, which connects 
these two. The evolution of the relative coordinate in 
each of these parts can then be described by two-by-two 
mappings of the relative coordinate and momentum. The 
composite mapping, built up of these building blocks, de- 
scribes the relevant dynamics. 



A. Mapping the free motion 

The trap potential V{x) (0) is flat in the central part 
of the trap and linear in the mirrors; in these regions the 
Hamiltonian (||) can be separated into CM and relative 
coordinates. 

Without the ions' interaction, the relative momentum 
p is constant in both the central part and in the mirrors 
and the distance x evolves freely as 



x{t) = x{0) + -t. 



(10) 



As the times spent in the mirror and in the central part, 
Tm and Tf — aTm, are fixed by the CM motion (see 
equations (||) and (||)), we define two mappings J- and M. 
that describe this free evolution of the relative coordinate 
in the central part and in the mirrors, respectively. They 
differ only in their time duration: 
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(11) 



^^"tA (12) 
P J 



Their interpretation is the following: if the ions have 
the relative coordinates {x,p), when the CM enters, e.g., 
the mirror, these will have evolved to {x\p') = A4{x,p), 
when the CM leaves this region again. 




FIG. 2: During the time Tk, when the two ions are in dif- 
ferent parts of the trap potential, they experience a different 
acceleration due to the mirror potential. Integrating the dif- 
ference over Tt results in a compressing change of the relative 
momentum (|l3|). 



B. The momentum kick approximation 

For that part of the motion when the two ions are 
in different parts of the trap potential, the CM and the 
relative motion are coupled for a time r^. This is the 
time between when the first ion crosses the kink between 
the flat part and one of the mirrors at ±^ and when the 
second one passes that point (see figure ||). During that 
time the ions are accelerated relative to each other by 
the trap potential. For the case of incoming ions, which 
is depicted in figure 0, this time is: 



Tk = 



P2/m 



After this time the momentum p2 will be the same, p'2 = 
P2, but p'l has become 



Pl+Tk 



dV 

dxi 



= Pi 



-F 



P2 



The two ions have nearly the same velocity, so ^ is ap- 
proximated by the velocity of the CM. 

If the energy of the two ions of about 4.2 keV differs by 
2 eV in the laboratory frame, their momenta will differ 
by about 1 a.u.: replacing ^ by then introduces an 
error on the order of 0(10^^). 

The relative momentum p — P'^^p^ now changes as (cf. 



P =P 



FM 
~2P~ 



P 



2fj. 



(13) 



This formula is valid for outgoing ions, too. Through 
the special geometry of the trap the momentum transfer 
and the time spent in the mirror are intimately related 
to each other. 

The relative distance x changes, too, but assuming an 
initial distance of, e.g., x — 10 a.u., only by a factor 
of 1 — 0(10"'^); this change of x will be neglected. We 
also assume that the CM momentum is not affected by 
this momentum transfer, which is true to the order of 
0(10-3). 

As the ions — in our approximation — do not move 
relative to each other during r^, the motion through the 
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FIG. 3: Sketch of the sequence of mappings during one period 
of the motion of two non interacting ions (|l|). The trap IS 
symbolized by the potential and the stylized pickup. 



kink region has the effect of an instantaneous momentum 
"kick" in the relative coordinate, induced by the coupling 
of the CM and the relative coordinate. This momentum 
"kick" always pushes the two ions together. We conse- 
quently describe the kink by the following mapping /C: 




(14) 




FIG. 4: Sketch of the trajectories of two non interacting ions 
on their way from the central field free region of the trap 
through the mirror back into the central part. The horizontal 
eads is time, the vertical axis denotes the position of the ions 
in the trap. Ion 2 is here the faster one. The initial relative 
momentum is p, after the first kick it is denoted by p and after 
the second kick by p' . The time when the ions' trajectories 
cross is denoted by i. 



with the same accuracy of 0(10 
persion (0). 



') as the linearized dis- 



C. The composite mapping — the dispersion of the 

trap 

With the mappings A4 and JC in hand we are now 
able to describe the motion of two non interacting ions 
in the trap. One complete period consists of a sequence 
of these mappings, cf. figure |^. We start the cycle in the 
middle of the central region, where in the experiment 
the pickup is located. The composite map, denoted by 
Vnon, then describes how the observed distance between 
the two ions changes between successive measurements. 
It has the following form, with JFi describing the free 

motion for a time ^ (cf. (|l] 



Vnon = ( 

2 



Q 



G 



^lC®Ti (15) 
).Fi(16) 



In the second line the effects of the mirror and the two 
adjacent momentum kicks are combined into the map- 
ping Q. We will later see that the main synchronization 
process can be understood from this part alone. The 
composite mapping Vnon finally simplifies to 



D. The mirrors' map: Q ^ K.® M®K. 

Now we will look at the above defined mapping Q (|1( 
From energy conservation and geometrical considerations 
we see that — without the ions' interaction — the total 
effect of the mirror is to turn around the direction of the 
two momenta pi and p2 , and therefore of p (see figure ^) : 



P'l 

P'2 




t Pl~P2 
P = 7\ = -P 



(18) 



The relative momentum changes by Ap — p' — p = —2p. 
Each of the two momentum kicks is proportional to the 



corresponding distance (g 



at the entrance and x' at 



the exit of the mirror; the intermediate and final mo- 
menta, p and p' , are therefore 



2p 

m. 



and 



P =P 



2/i 

m. 



(19) 



From the above equations (|T^) and (|l9|) we derive a rela- 
tion between the distances x and x' and the momentum 



p = -p 



Vn 




2T„ 



(1 - a)p 



P 



(17) 



From energy conservation we know that the momenta pi 
and p2, and consequently p, have to be the same, when 
each of the independent ions has completed the period. 
The increase in distance, already calculated from the dis- 
persion ^ (0) is also reproduced by Vnon- 

Our mapping model therefore describes the evolution 
of the relative distance between two independent ions 



X + X = p 

M 

This relation can also be obtained by evaluating Q = 
K® M®1C directly: 




-P 



(20) 



It should be noted that with the experimental param- 
eters of table |, the change of x in the mirror is about 
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a factor of 10 bigger than Ax after one complete period 
(^. Most of the mirror's effect, however, is compensated 
for by the motion in the central part. At = 0, which 
corresponds to a = 1, both contributions exactly cancel 
each other. 

In figure ^ the time i denotes the crossing of the two 
ions' trajectories. It is defined by 



and evaluates to 



i = i{x,p) 



2fIX - TmP 



(21) 



The limits for 



ztoo 01 p — are t — For 



finite p and - < 0, < t < for - > 2k we have 

^ p ' — — 2' p — fj, 

2k- < i < Tra- In the intermediate region t is either 
negative or bigger than Tm, with a singularity at 
In that case there is no collision inside the mirror, as 
( pT| ) is only defined in the mirror region. 

Equation (^ ) is meaningful only as long as both ions 
spend time together in the mirror: if the first ion has 
already left the mirror when the second one arrives, the 
above treatment is not valid. This happens when the ions 
are further away from each other than 



P L 
M a 



(22) 



Without interaction the two ions' trajectories cross twice 
during a given period. This happens either in the mirror 
or in the central region. The probability wp for the latter 
is given by the ratio between the relative distance for 
which there is no solution of ( pl| ) inside the mirror, ^i^, 
and the distance Xmax (^2|), for which equation (^) is 
defined. With (||) it evaluates to 



Mp _ Ap 
~iIP ~ ^ 



(23) 



As p ^ P, a collision of the two ions in the central part 
is a highly unlikely event. With the parameters given, 
wp is on the order of only 0(10"'^). 

The different possibilities can be visualized by shifting 
the two (parabolic) trajectories of the ions in figure ^ 
against each other: as they are plotted in figure ^ we 
have < i < 3k _ This is case (i) of figure ||. If now the 
trajectory of ion 1 is moved to the right, the crossing time 
t will slide to the left, until both trajectories intersect 
at their entrance into the mirror, which means t = 0, 
depicted in case (ii) of figure ||a. When trajectory 1 is 
shifted further to the right, case (iii), the ions do not 
cross inside the mirror (cf. the explanation above and 
equation (p3|)). The crossing point reappears ai i — T,„ 
(iv), i.e., when the ions leave the mirror, and from there it 
proceeds back in the direction of 2k (y). These processes 
are the framework into which we later incorporate the 
ions' interaction. 



Fig. ||(b) plots the relative momentum p{t) during the 
course of Q. In all cases the final momentum is p' = —p. 
For cases (i) and (ii) the momentum changes sign only 
with the second kick, while for cases (iv) and (v) it is 
reversed already with the first kick. When the crossing 
takes with i < 2k the second kick is the larger one, while 
ii i > 2k ^ the first kick is stronger. In all these cases the 
two kicks have opposite direction. 

In the above cases x(t) changes sign in the mirror at 
t = t. In the special case (iii) x does not change its 
sign and both kicks work in the same direction; here the 
momentum in the mirror is less than before and after it. 
In this regime the faster ion enters the mirror first and 
leaves it second. 

It should be emphasized that, without interaction, the 
relative momentum in the central part of the trap is the 
same after each period. This is the important conserved 
quantity. The relative momentum in the mirror, on the 
other hand, is a measure of the time delay or phase lag 
of the two ions on their orbit. 



IV. ADDING THE ION-ION INTERACTION 

We now incorporate the (repulsive) interaction be- 
tween the two ions into our model. We will see that 
the repulsive interaction can be modelled as a time delay 
and we will then introduce this delay into the mappings. 

A. The time delay 

When both ions are in one of the mirrors or in the 
central part, the relative coordinate is decoupled from 
the CM motion (see section III A). The free motion in 
the relative coordinate is now modified by the interaction 
W{x) (p^). Any collision of the ions is elastic due to the 
conservative Coulomb interaction. As we are only inter- 
ested in the final values of x and p, but not in the actual 
solution of the equations of motion, we adopt a central 
idea of scattering theory: the whole effect of the ions' col- 
lision is described by a phase shift, or, for our treatment, 
a time delay r This delay modifies the propagation 
time for the freely evolving relative coordinate: 



x{T) ^ x{0) + ^{T ~ t) 



(24) 



In this ansatz T stands for either Tm or Tf. The time 
delay t is clearly a function of the interaction potential 
and the initial distance and momentum. It can be posi- 
tive or negative; if the relative motion is delayed we have 
T > 0. 



B. Modifying the mappings 

With the interaction added we now use (|2j) instead of 
( pO| ) to map the relative coordinate in the mirrors or the 



6 




trap's central part, respectively. 

If the initial and final distances are not in the asymp- 
totic region, where W{x) — 0, we have to adjust p' , too, 
to conserve the energy of relative motion. But as our 
treatment does not depend on a special form of W, we 
will for the following assume that the collision between 
the two ions takes much less time than Tm or Tf and that 
before and after the collision the interaction between the 
ions can be neglected. 

A time delay t„i in the mirror region may be incorpo- 
rated directly into (p^). It modifies the relative coordi- 
nate: 



M' 




i.p 



P 



(25) 



The mapping now is a functional of the still unspecified 
delay r™ = Tm{W;x,p), which itself depends on the ac- 
tual form of the interaction potential. 

In the non interacting case the mapping in the central 
part is described by Recall that the ratio of Tj 

and Tm defines the parameter a, which in turn is related 
to the dispersion of the trap (^. It will be convenient, 
therefore, to incorporate a time delay Tf into a modified 
a', defined by: 



- Tf {a - ■^)T,n =: a'T^ 
The mapping T is consequently modified to: 




i'T„ 



P 



(26) 



(27) 



During the short interval of the momentum kick we as 
sum e the ions do not move relative to each other (see sec- 
tion 



[II B). Also the trap potential is much stronger than 
the weak ion-ion interaction. Consequently the mapping 
K ( p^ ) remains unchanged. 



With these modifications the composite mapping for 
the whole period, starting and ending at the pickup in 
the center of the trap, becomes a functional of the delays 
and Tf (cf. (|l|)): 

V ^T[®K®M' (E)1C®T'®K®M' ®K®T[ (28) 

2 2 

Here T't denotes the propagation through only half of 

2 

the central part, i.e., between the pickup and one of the 
mirrors: T' — T\ ® T'x ■ 



V. LINEAR STABILITY ANALYSIS 

Now that we have set up all parts of the model we will 
quantify the stability of the composite mapping V' ( |2^ ) 
and determine the conditions under which synchroniza- 
tion of the two ions' motion occurs. 

We follow the usual lines: the mappings are linearized 
around the initial values {x,p) by calculating the Jacobi 
matrix of partial derivatives: 



dx' 


dp' 


dx 


dx 


dx' 


dp' 


dp 


dp 



J = J{Tm,Tf) 



If the eigenvalues A of J are complex with a length of 
I, then the mapping is stable and the ions' motion syn- 
chronized; if the A are real it is unstable and the motion 
diffusive g. According to equations (|l|), (|2|) and (|2^) 
the eigenvalues are functions of the trap parameters Tm, 
a and and of the time delays t„i and Tf. 



A. The elementary maps 

Through the linearization the time delays become con- 
stants. With the constant and Tf, and, consequently. 
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the constant a', all the maps involved are linear and their 
Jacobians can be represented by 2 x 2 matrices. They 
have the following forms: 



1 

1 

1 1^ 







1 



.2lL 1 



(29) 
(30) 
(31) 



The determinant of each of these matrices is equal to 1. 
Moreover, they all have a double eigenvalue of 1 and a 
fixpoint at {x,p) = 0, independent of the values of T^, 
/i, a' and r^. 

The composite mappings are represented by the cor- 
responding products of the elementary matrices: G = 
KMK and P = F i KMKFKMKF i , where Fi corre- 

2_ . 2 2 

sponds to T'l (cf. (^). Since F, M and K all have a 

fixed point at {x,p) = 0, it follows that G and F must 
as well. 



B. Stability of the complete mapping V 

The Jacobian for the complete period can be grouped 
into the modified mirror G = KMK and the field free 
central region F: P = FiGFGFi. The stability cri- 
terion can, due to this repetitive structure, already be 
inferred from Pi = Fj.GFi. Its explicit form is. 



-l + 2e(l 



^(1 



a){l- 
f 2e(l 



a) 



where we have defined e = ^ . 



(32) 



The determinant of P i is 1. In section HID we showed 



that the ions collide either in the mirror or in the cen- 
tral part; therefore one of the delays vanishes. These two 
cases correspond to either a' = a or = 0, respectively. 
For a' = a and Tm = the non interacting case is re- 
covered: P = PiPi = Pnon, where Pnon IS the matrix 

2 2 ' ^ ^ 

corresponding to the mapping of equation (^11 



1. Collisions in the central part 

First we will consider the case of = 0, i.e., the 
collisions take place in the central part. As discussed 
above this case is rare. Then P = Pj. x Pi reduces 
to the dispersion of the non interacting case (p7|), but 
with the externally defined a replaced by the modified 
a' = a~ Y~- Consequently Pi has a double eigenvalue 
of —1, i.e., the distance between the ions grows linearly 
in time according to equation (R|) with a replaced by a' . 



Which ion is faster is determined by whether a' is bigger 
or smaller than 1. 

As in the non interacting case, the eigenvalues of the 
composite map for the whole trap are A1/2 = 1: the ions 
separate linearly in time, if the modified dispersion is not 
exactly a' — 1. 

For a < 1, which is the experimentally observed con- 
dition for bunching, a negative delay Tf may bring the 
trap into the regime where the "effective" dispersion a' 
becomes 1 and both ions have the same period. Colli- 
sions in the central part may therefore synchronize the 
ions by modifying the dispersion of the trap: their dif- 
fering momenta remain unchanged, but the distance is, 
due to the vanishing dispersion, the same every time they 
pass the trap's center. The possibility of this type of syn- 
chronization is small: first collisions in the central part 
are rare events (see equation (^3|)) and second, with the 
eigenvalue of 1 the synchronization, if it occurs, is not 
stable against perturbations. 



2. Collisions in the mirror 

In the other, much more frequent case, where the ions 
collide in the mirrors, i.e., a' = a and Tm 7^ 0, the eigen- 
values have the form 



Ai/2 = -1 + 2e(l - a) ± 2^e'^{l - af - e(l - a) . (33) 

The unmodified a describes the trap's original dispersion. 
It is now convenient to define another parameter 



7 = e(l — a). 
Then the eigenvalues (^3|) simplify to 



Ai/2 = -1 + 27 ±2^7^ 



7, 



(34) 



(35) 



The argument of the square root, 7^ — 7, is negative for 
< 7 < 1. In that region the eigenvalues are conse- 
quently complex with a length of 1 and they are real for 
7 < 0. 

The important parameter, as we see, is composed of 
the time delay Tm in the mirror and the dispersion of 
the trap, described by a (0). The sign of this parameter 
7 determines if the relative motion of the two ions is 
bounded or not: for 7 > 0, which requires either t„j > 
and a < 1 or Tm < and a > 1, the eigenvalues A are 
complex and the ions' motion is synchronized. For 7 < 
the eigenvalues are real and the ions separate faster than 
without interaction. This happens for t,„ > and a > 1 
or for Tm < and a <1. 

For very small 7, i.e., ^ T,„ or a w 1, we can 



neglect 27 compared to 1 in (35) and approximate the 
square root by ^77^ — 7 ~ V~7- The eigenvalues are 
then of the form 



lim Ai /2 = — 1 ± 

7->0 ' 



(36) 
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r„ > 


r,„ < 


> Tf 

— r\, ^ 


7 > 
A complex 
synchronization 


7 < 
A real 
difltusion 


< 

Q > 1 


7 < 
A real 
difltusion 


7 > 
A complex 
synchronization 



(a) case (i) ^ 

kXiit ) 



TABLE II: Overview of the eigenvalues A of the composite 
mapping matrix Pj_ for collisions in the mirror and the 
corresponding behavior of the ions. For 7 = 0, i.e., Tm, = 
or Q! = 1, which corresponds to either collisions in the central 
part or no interaction at all, the eigenvalues are —1 and the 
ions' separation in phase space grows linearly in time. 



Consequently the transition between real and complex 
eigenvalues is very abrupt and the stability changes com- 
pletely as the delay or the dispersion changes its sign. 
Table gives an overview of the behavior of the eigen- 
values for the different values that a and t„j can take for 
collisions in the mirror. 

The stability analysis is consistent with the experimen- 
tal observation and the findings of reference |^ that syn- 
chronization occurs in the region of ^ > 0, which in 
our parametrization corresponds to a < 1, or accord- 
ing to to 2Eap > (|^): to achieve bunching, ions 
with higher energy must have a longer period. But this 
criterion alone is not enough; it has to be accompanied 
by a time delay in the mirror, which in turn requires a 
repulsive ion-ion interaction. 

For the other scenario, where the mapping is stable 
for a > 1, the time delay has to be negative. This is 
possible if the ions have an attractive potential, e.g., if 
their charges have opposite signs, but then the electro- 
static field of the mirrors could keep only one of the ions 
trapped; their motion would not even have a chance to 
synchronize. For two or more equal ions only a repul- 
sive Coulomb interaction is possible; any other interac- 
tion like the attractive van der Waals interaction is much 
weaker and can not compete with the Coulomb repulsion. 
A negative time delay is also possible when the ions' re- 
pulsion becomes comparable to the guiding mirror field: 
then a repulsion between the ions can lead to a nega- 
tive delay — the ions "bounce off" each other. But in 
this case the ions' repulsion would dominate the dynam- 
ics and it is questionable if the very weak trap potential 
would be able to actually trap the ions. 



VI. EXPLAINING THE SYNCHRONIZATION 
MECHANISM 

From the stability analysis we identify two different 
mechanisms that can synchronize the motion of two ions: 
when the ions collide in the mirror, the relation between 
the momentum kicks is unbalanced, and when collisions 
take place in the central part, they modify the effective 




FIG. 6: (a) The ions' positions in the mirror vs. time for 
case (i), where t < Tm/2 and for case (v), where t > Tm/'i 
(cf. figure (b) The corresponding evolution of the relative 
momentum. A time delay in the mirror leaves the ions 
less time to separate after their collision. Consequently the 
second momentum kick is smaller than without interaction 
(gray lines). 



dispersion. We will now have a closer look at the micro- 
scopic dynamics in these two regimes. 



A. Collisions in the mirror 

To explain how the collisions in the mirror synchronize 
the ions, it is sufficient to look at the modified mirror 
map Q' = K, ® M' ® K,. Without the interaction, Q 
exactly turns around the relative momentum of the two 
independent ions: the two momentum kicks add up to 
— 2p. During the interval in the mirror the ions approach 
and reseparate. When the repulsive interaction is added 
the ions experience a time delay and, consequently, have 
less time to separate: the second momentum kick ( [T^ ) 
will be smaller than without interaction. This is depicted 
in figure^. The two sets of trajectories correspond to the 
configurations (i) and (v) of figure ^ which shows the 
non interacting case. Now the two kicks, which both are 
proportional to the corresponding relative distance, will 
not add up to —2p any more and the relative momentum 
after the second kick will be either smaller (figure ^ case 
(i)) or bigger (case (v)) than — p, depending on whether 
the slower or the faster ion enters the mirror first. In 
either case momentum is effectively transferred from the 
second onto the first ion, due to the modified coupling 
between the relative and the CM motion. 

If the first ion is the faster one it will gain additional 
energy through this momentum transfer. Because the pe- 
riod increases with energy, the distance between the two 
ions will shrink until the slower one overtakes the faster 
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one. Then momentum is transferred from the faster ion, 
which is now the second one, onto the slower first one, un- 
til the first ion again will be the one with the higher mo- 
mentum. This process repeats itself over and over, with 
the ions effectively orbiting around each other in their 
relative coordinate, stably interlocking and synchroniz- 
ing their motion. 

In the non interacting case the momentum kicks and 
the time spent in the mirror are intimately related (see 
(pT|)). This balance is now offset by the ion-ion inter- 
action. In this way the very small repulsion is amplified 
to counterbalance the trap's dispersion and suppress the 
separation of the ions. 

B. Collisions in the central part 



0.6 

13 



-0.6 



-4 -2 2 4 -4 -2 2 4 
x[10^ a.u.] 

FIG. 7: Poincare section for the composite mapping V for a 
constant Tm = 3 x 10* a.u. Panel (a) shows the bunching and 
(b) the diffusive behavior. The arrows denote the direction of 
the phase space flow. 




When the ions collide in the central part of the trap the 
mech anism is completely different. As shown in section 
III D , in this regime the faster ion enters the mirror first 
and leaves second (see figure ||(iii)). For this mode to 
be stationary the dynamics has to be symmetric in both 
mirrors; consequently the ions have to exchange their 
roles in the collision in the central part of the trap: the 
first slower ion is accelerated by the faster second, which 
comes from behind. With this mechanism, which was 
dubbed "bricking motion" by Zajfman etal 1^, each of 
the ions is the faster for one half period and the slower 
one on the other side of the trap. This way both ions 
have the same energy, and therefore the same period on 
average, even for a ^ 1. 



VII. ILLUSTRATIONS AND EXTENSIONS 

In the previous sections we have derived the conditions 
for stability and explained the underlying mechanism of 
our mapping model. Now we will check the validity of our 
description with regard to the two central simplifications 
— the instantaneous ion-ion interaction and the constant 
slope of the mirror potentials. 

A. Poincare maps for constant Tm 



Poincare section becomes hyperbolic, see panel (b): the 
ions separate with a speed that increases with their dis- 
tance. 

Note that the phase space flow in figure |^(a) is in the 
direction opposite to that of an isolated particle in an 
external harmonic potential. 

For a constant Tm and a < 1 , using the linearized map- 
ping Pi ( p2| ) we can derive a conversion factor between 
the maximal distance x^ax and the maximal momentum 
Pmax = ^7^°°' of a Poincare ellipse. It evaluates to 



2fj. 
T 

2^ 



1 — a 



{l-af 



7-7 



(37) 
(38) 



The parameter 7 (|3J) is the same as defined in the con- 
text o f the stability of the bunching mode (cf. section 
VB2), but the sign of the argument of the square root 
is reversed. For the Poincare sections shown above with 
= 3 X 10^ a.u. we get /? = 7.6 x 10^ a.u. 
From equation (|3^) we see that f3 decreases when the 
trap is brought closer to the transition at a = 1: the 
maximal momentum increases until the point that the 
ellipses "break apart" and become hyperbolas. 



In section VIA we explained that in the bunching 
mode the ions oscillate around each other. This "or- 
biting" of the ions can be conveniently depicted in a 
Poincare section of the repeated mapping V' (28) by plot- 
ting the relative coordinate and momentum after each 
iteration. 

In section^ we linearized 7^', which resulted in a con- 
stant Tm ■ Figure shows the results of the repeated map- 
ping for an arbitrarily chosen = 3 x 10^ a.u.: Panel 
(a) is calculated for a = 0.956 based on the experimental 
values given in table ^, resulting in ellipses. 

If we increase a to 1.046 = 1/0.956 into the unsta- 
ble regime, with all the other parameters unchanged, the 



B. Trajectory calculations: non instantaneous 
delay 

With the simple mapping model we were able to derive 
the condition for a stable bunch and to give a microscopic 
explanation for the coupled motion of the two ions. In 
order to verify that the discrete description of the rela- 
tive coordinate implied by the mapping is accurate, we 
compare our results to trajectory calculations, i.e., to the 
exact solution of the equations of motion of the two ions 
from the full Hamiltonian (||). 

To allow the ions to pass by each other in our one di- 
mensional model we replace the Coulomb repulsion jij 
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0.25 0.5 0.75 1 
t/T 

m 

FIG. 8: Solution of the equations of motion for two ions 
for one period T with and without interaction. These plots 
should be compared to the results of the mapping (figure |^ . 
Panel (a) plots the positions: the distance of the two ions 
( ) from the CM's trajectory ( ) is multiplied by a fac- 
tor of 300. The horizontal lines at a;i/2 = ±1.89 x 10^ a.u. 
denote the beginning of the mirror regions, (b) shows the rel- 
ative momentum between the ions with ( ) and without 

( ) the interaction. The temporary slowing down due to 

the interaction is clearly visible. 




-4 -2 2 4 -4 -2 2 4 
X [10^ a.u.] 

FIG. 9: Poincare sections from trajectory calculations for the 
stable regime of q = 0.956 (a) and for the unstable regime 
with a — 1.046 (b). The lower row is the result of the repeated 
mapping V with a non constant Tm from (^) for the same 
values of a as (a) and (b): (c) shows the stable and (d) the 
unstable regime. 



between the ions by a so called " softcore Coulomb inter- 
action" [ITol: 



W{x) = 



1 



The "softcore parameter" or "impact parameter" c? is a 
measure of how close the two ions have to come when 
passing each other; it describes the diameter of the ion 
beam in the trap potential. For the trajectory calculation 
we set this value to d = 10 /im. This value is much 
smaller than the real beam, but still the ions' interaction 
is weaker than the trap potential by about four orders of 
magnitude. 

Figure ^ we show the result of a trajectory calculation 
for one period T for the initial values x — —10^ a.u. 
and p — 1 a.u. All other parameters were set to the ex- 
periment's values of table||, resulting in a = 0.956. Panel 
(a) shows the CM's position vs. time as a broken curve 
and the two ions' coordinates. Their distance from the 
CM has been increased by a factor of 300 for clarifica- 
tion. Case(i) of figure ||(a) is reproduced during the time 
interval O.IT < i < OAT and, symmetrically in the other 
mirror, during 0.6T < t < 0.9T. Panel (b) compares the 
relative momentum with (solid line) and without (broken 
line) interaction: the "discontinuities" of the momentum 
kicks are clearly seen, cf. figure ^. Without interaction 
the momentum is reversed after each mirror and returns 
to its initial value: p" = —p' = p after one full period. 



The main difference between the trajectory calcula- 
tions and the mapping description is that in our map- 
ping approach the ions' collisions are instantaneous and 
confined to the mirror. In the trajectory calculations the 
ion-ion interaction takes place anywhere in the trap on a 
time scale comparable to Tm for small relative momenta. 
This can be seen in figure |8|(b) where the relative mo- 
mentum deviates visibly from the non interacting case 
for about half of the time spent in the mirror. 

The trajectory shown in figure ^ confirms that the pro- 
cesses of figures H and |[ which we used to explain the 
mechanism, exist. In order to confirm both the bunching 
condition and the dominance of the collisions in the mir- 
ror, we created Poincare sections from the trajectory cal- 
culation, too. The stability of the synchronization shows 
up in the overall elliptic or hyperbolic structure of the 
orbits. If the ions are synchronized dominantly through 
delaying collisions in the mirror, as we propose, the or- 
bits will be (deformed) ellipses, whereas the alternating 
"bricking motion" (cf. section VI B ) will lead to a dis- 
crete structure. 

The Poincare sections from the trajectory calculation 
are shown in figure ^ in panels (a) and (b) . In the sta- 
ble regime, i.e., for a < 1, shown in (a), the orbits are 
deformed ellipses: the ions orbit around each other. The 
strength of the interaction, and therefore the local cur- 
vature of the orbit, varies with the relative distance and 
the momentum, but it can still be summarized by a time 
delay. 
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Only for small distances and momenta does the struc- 
ture of the Poincare section become more complicated: 
the fixpoint at = is surrounded by a chaotic 

layer. When the ions are very close together the effective 
delay changes so much between successive collisions that 
the quasicontinuous orbits are broken up and consecutive 
points are distributed randomly in the central region of 
the Poincare section. 

The agreement between the Poincare sections from the 
mapping and from the trajectory calculations can be im- 
proved, when the constant of figure ^ is replaced 
by a momentum dependent delay. That the delay de- 
pends only on the relative momentum is motivated by our 
ansatz that the ions reach the asymptotic region again af- 
ter the collisions. Then the delay does not depend upon 
the initial and final distance, but only on the speed with 
which the ions pass each other. In figure |^(c) we show 
the result of the mapping with the ansatz 



c 



1 



(39) 



The functional form and the parameters — C ~ 3 x 10^ 
a.u. and po = 1 a.u. — were chosen to give orbits of a 
comparable structure and size as the trajectory calcula- 
tions of panel (a). 

In figure ^(b) we increased a into the unstable regime 
to a = 1.046. Then the trajectory calculation becomes 
unstable, just as the mapping, except for a small island 
around {x,p) = 0. Even this island vanishes when a is 
increased further. The stability condition is a < 1 as in 
the mapping, but the transition is less abrupt than with 
the mapping approach. 

Panel (d) plots the corresponding Poincare section of 
the mapping with the momentum dependent delay of 
equation ( p9[ ) . Again, as was the case with panels (a) and 
(c) in the bunching regime, the overall structure is com- 
parable to the trajectory calculation of (b). The small 
stable island does not exist, though. 

The Poincare surfaces of section provide a further in- 
dication that the mapping ansatz is an excellent repre- 
sentation of the full trajectory dynamics with respect to 
both the dominant synchronization mechanism and the 
stability condition. 



C. Non— linear mirror fields 

Above, we modelled the trap with a constant gradient 
mirror field (|^). This allowed us to decouple the relative 
coordinate from the CM motion. We now have to ensure 
that our results and explanations are also valid when the 
mirror field is not linear in the experiment. 

Two conditions must be fulfill ed for sync hroniz ation to 
occur in our model (cf. sections VB 2 and VIA): (i) the 



total dispersion of the trap has to be positive and 
(ii) the interaction has to delay the ions in the mirror, so 
that the relative momentum is not exactly reversed. 



It is easy to see that the overall positive dispersion is 
not restricted to constant mirror field configurations; any 
mirror field that increases slower than harmonic will do. 
The mirror only has to be "long" enough so that it can 
compensate the negative dispersion of the field free part 
of the trap. 

If the mirror potential is not linear it nevertheless re- 
flects the individual momenta of two non interacting ions. 
The relative and the CM motion can not be separated 
any more, but the evolution of the relative coordinate 
— and consequently the momentum reflection — can be 
described by an alternation of infinitesimally small kicks 
and free evolutions. 

The ions' interaction is then added into this succession 
of mappings as a set of infinitesimal delays, each modify- 
ing the balance between the adjacent kicks. The overall 
effect is the same as with the linear slope: when the ions 
repel each other they leave the mirror with a smaller dis- 
tance and the second ion transfers energy onto the first. 

Consequently the conditions for synchronization and 
the explanations given remain unchanged for an arbitrary 
mirror field. 



D. Connection to previous explanations 

The approach presented here has several aspects in 
common with the one recently given by Strasser etal . 
In fact the propagation matrices for the non interacting 
ions are identical. In the Strasser model the composite 
mapping (^7|) derived from these matrices is cast into an 
effective reduced mass m* (see equations (6) and (8) of 
reference |^). When the ions repel each other, m* has 
to be negative for the relative motion to be bound. This 
leads to the bunching condition of a positive dispersion 
of the trap, i.e., > 0. Both the trap's dispersion and 
the ions' repulsion are described in a mean field ansatz; 
the Strasser model is therefore very robust with respect 
to the exact form of the trap potential or the details of 
the ions' interaction, but can not give an explanation of 
the underlying microscopic dynamics. 

Both the Strasser and our model require a < 1 for 
bunching, but whereas in our approach any delay r™ > 
is sufficient to couple the ions, the Strasser model requires 
a minimum density in the bunch, i.e., a minimum inter- 
action strength, because it confines the test ion to stay 
inside the bunch. 

By comparing equation (10) of reference with our 
equation ( p2[ ) we can relate the static bunch's parameters 
to an effective time delay fm- 



(1 



(40) 



The bunch consists of N ions of charge q and mass m 
and has a radius Rq; eg is the dielectric constant. 

From the phase space flow of the Poincare sections we 
see that in the diffusive regime the faster ions pass the 
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pickup in the front of the bunch. This was already ob- 
served in In the bunching regime, though, this or- 
dering is not simply inverted. The Poincare sections are 
symmetric with respect to the coordinate and the mo- 
mentum axis, consequently each of the ions is on average 
in front for half of the time and has the higher momentum 
for half of the time, when the CM passes the pickup. 

VIII. SUMMARY AND CONCLUSIONS 

In this paper we have set up a model to describe the 
motion of two ions in an "ion trap resonator" , an ion trap, 
which essentially consists of two spatially separated elec- 
trostatic mirrors. Our one dimensional model describes 
the evolution of the ions' relative coordinate and mo- 
mentum with three simple mappings: one, each, for the 
free evolution in the central field free region of the trap 
and in the mirrors, while the third mapping describes the 
connecting kink between these regions: here the relative 
coordinate is coupled to the center of mass motion. The 
coupling gives rise to an instantaneous compressing kick 
in the relative momentum. The interaction between the 
two ions is summarized by a time delay, which shortens 
or lengthens the time for the free evolution either in the 
mirror or in the central part. 

The three mappings allow us to perform a linear sta- 
bility analysis of the composite mapping for the whole 
trap. Without delay, i.e., without interaction, the map- 
ping model reproduces the behavior of two independent 
ions; with the delay added we identify the criteria for 
the stable synchronized motion: the trap has to be oper- 
ated in a regime where the period of the ions in the trap 
increases with their energy and the ions have to be de- 
layed by their repulsive interaction, when they cross their 
paths in the mirror. This confirms the experimental find- 
ings that the dispersion of the trap has to be positive. 

Based on these simple building blocks we are able to 



describe, how the interaction modifies and couples the 
motion of the two ions in the bunching regime: without 
the interaction ions with different energies separate lin- 
early in time due to their different periods in the trap. 
A time delay in the mirror, which describes the repulsive 
interaction, now modifies the balance of the coupling be- 
tween the relative and the CM motion at the kink, so 
that during each pass through the mirror momentum is 
transferred from the second onto the first ion. This ad- 
ditional energy increases the first ion's period and lets it 
fall back against the other one until they have exchanged 
their places. Now the energy is shuffled back onto the 
other ion, until they exchange their position again and 
the circle is completed: in the stable region the ions os- 
cillate around each other, constantly transferring energy 
back and forth between them. 

Trajectory calculations confirmed that the ions are 
synchronized dominantly by the microscopic process of 
collisions in the mirror that we had derived from the 
mappings together with an instantaneous interaction. 
This continuous description reproduces the stability con- 
ditions, too. We also verified that our model is valid for 
arbitrary forms of the mirror potential. 

Our model of two ions can also be applied to the dy- 
namics of a test ion inside a bunch. The microscopic 
picture of the ions' motion, which was developed in this 
paper, then allows us to identify the important regions of 
the trap and to explain how they shape the bunch. This 
will be the subject of a forthcoming publication 
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